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ABSTRACT 

Some of the difficulties in determining the underlying physical properties that are relevant for 
observed anomalies in microlensing light curves, such as the mass and separation of extra- 
solar planets orbiting the lens star, or the relative source-lens parallax, are already anchored in 
factors that limit the amount of information available from ordinary microlensing events and 
in the way these are being parametrized. Moreover, a real-time detection of deviations from an 
ordinary light curve while these are still in progress can only be done against a known model 
of the latter, and such is also required for properly prioritizing ongoing events for monitor- 
ing in order to maximize scientific returns. Despite the fact that ordinary microlensing light 
curves are described by an analytic function that only involves a handful of parameters, mod- 
elling these is far less trivial than one might be tempted to think. A well-known degeneracy 
for small impacts, and another one for the initial rise of an event, makes an interprediction of 
different phases impossible, while in order to determine a complete set of model parameters, 
the fundamental characteristics of all these phases need to be properly assessed. While it is 
found that the wing of the light curve provides valuable information about the time-scale that 
absorbs the physical properties, the peak flux of the event can be meaningfully predicted only 
after about a third of the total magnification has been reached. Parametrizations based on ob- 
servable features not only ease modelling by bringing the covariance matrix close to diagonal 
form, but also allow good predictions of the measured flux without the need to determine all 
parameters accurately. Campaigns intending to infer planet populations from observed mi- 
crolensing events need to invest some fraction of the available time into acquiring data that 
allows to properly determine the magnification function. 
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1 INTRODUCTION 

An efficient detection of planets by means of gravitational mi- 
crolensing requires sufficiently accurate predictions of the under- 
lying ordinary light curve against which the pl anetary deviations 
need to be identified jPominik et alj2007l . l2008l) . An optimal mon- 
itoring strategy in order to maximize the scientific return moreover 
profits strongly from the proper determination of a full set of model 
parameters as early as possible, and the latter becomes a require- 
ment for finally assessing the planet detection efficiency and draw- 
ing conclusions about the planet population. 

It is kno wn that early - stage event prediction suffers from 
degeneracies jAlbrow|[2004h . in particular the peak flux is hard 
to assess, while the event time-scale tE, required to relate the 
observations to the underlying physical properties, can remain 
strongly uncertain even after the event has been observed over its 
full course, but not covered well enough jWozniak & Paczvnskil 
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Il997h . The lack of determinacy of Ie is particularly apparent for 
strongly-blended events, which in fact comprise the full sample 
for observations to wards neighbouring galaxies such as M31 (e.g. 
lBailIonetalJl993h . 

Understanding of ordinary light curves (comprising single 
point-like source and lens stars), their optimal parametrization, and 
the apparent degeneracies and ambiguities, is also crucial and use- 
ful for modelling events that involve anomalies, and drawing con- 
clusions about e.g. stellar binaries, stellar masses derived from par- 
allax measurements, and planets. Sparse event coverage in criti- 
cal regions is even prone to lead to ordinary events allowing for 
multiple minima of the (least-squares) hypersurface, this fact 
being further complicated by the presence of potential outliers in 
the data and th e application of robus t-estimation techniques to deal 
with these (e.g.l Pominik et aOlOOTh . 

In this paper, different phases of ordinary microlensing events 
along with their characteristics are identified, and it is shown how 
feature-oriented parametrizations can be used to quantify the be- 
haviour while avoiding parameter correlations. Moreover, some 
fundamental requirements for a monitoring strategy that allows 
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to meet the goal of determining either the observed flux or the 
corresponding magnification are discussed. The basis for this is 
formed by approximate relations between the magnification and 
the angular separation between lens and source that make the 
light curve independent of the impact parameter. Together with 
some arising degeneracies, t hese have already been identified by 
IWozniak & Paczvnskil ( Il997h . but here some new light is shed on 
the implications and the focus is on different aspects that are emer- 
gent right now, which leads to arriving at some new and different 
conclusions. 

While Sect. [2] reviews the general properties of ordinary mi- 
crolensing events and the canonical parametrization, and Sect. [3] 
presents the advantages of an alternative parametrization that is 
oriented towards the observable characteristic features, Sect.|4]dis- 
cusses the predictability of ordinary microlensing events or the lack 
of it. A summary and final conclusions are provided in Sect.|5] 



2 THE CANONICAL TREATMENT OF ORDINARY 
EVENTS 

A microlensing event results if two stars happen to be closely 
aligned as seen from Earth, where this alignment is quantified 
by the unique characteristic scale of gravitation al microlensing, 
namely the angular Einstein radius ( lEinsteinlll936[) 



4GM TTLS 



1 AU 



(1) 



where M denotes the mass of the foreground lens star, G is the 
universal gravitational constant, c is the vacuum speed of light, and 
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stands for the relative source-lens parallax, with Dl and Ds be- 
ing the distance from Earth to the foreground (lens) star and the 
observed background source star, respectively. With u 9e denoting 
the angular separation between lens and source star, gravitational 
bending of light by the lens star yields an observable magnification 
of the source star as an analytic function of u, which reads 



A{u) = 



(3) 



Stellar kinematics implies a non- vanishing relative proper mo- 
tion fi between lens and source star, so that u is a function of 
time and the magnification A(u) describes a characteristic light 
curve. While in the earlier history of gravitational microlensing 
several different even t time-scales have been used (e.g. iGrieslil99ll : 
Ide Ruiula et al.ll9^lh , Ie = Se/|mI has emerged as a popular con- 
venient choice. 

For uniform proper motion, i.e. constant fi, the dimensionless 
separation u can be expressed by means of three quantities that 
form the parameter vector p = (uq, to, ^e), so that 



u{t;Uo,to,tE) = -l/Mo + 



(4) 



where mo 9e is the smallest angular separation, encountered at 
epoch to, while the source moves by an angul ar Einstein radiu s 
relative to the lens within the time-scale tE ("Paczvnskil 1 1986h . 
Whereas uq and to only define the position of the trajectory of the 
source relative to the lens, and therefore do not carry any informa- 
tion about the relevant physical properties that determine the mag- 
nification, which are the relative parallax ttls, the relative proper 
motion /i, and the lens mass M, all these are convolved into the 



event time-scale ts- It is in fact its relation to the physical event 
properties that makes the parameter tE a preferred choice amongst 
possible time-scales. As illustrated in Fig.[T] ordinary light curves: 
(1) are symmetric with respect to a peak at epoch to, (2) reach a 
peak flux there, (3) approach a baseline flux for times far away 
from the peak, and (4) show characteristic inflection points. How- 
ever, the parameter to is the only one that directly relates to the 
characteristic features of the light curve, while all the others are not 
a proper reflection. This is not at all favourable for modelling, and 
instead being able to essentially read off the model parameters from 
the collected data would ease life a lot. 

The phenomenon of gravitational microlensing leads to a char- 
acteristic magnification A(t;p) as a function of time t and the 
model parameters p that describe the lens-observer-source geome- 
try, the lens properties, and the source brightness profile. The ob- 
served flux f'*' (t; p) for a given site and passband — denoted by 
the multi-index k — , however, furthermore is a linear function of 

(k) 

the intrinsic flux of the observed source star Fg , which is magni- 
fied, and a background flux F^''\ where (e.g. lAlbrow et alfeOOOh 



F(^)(t;p) = Ff^A(t;p) + 4*' 



(5) 



This allows us to isolate these two parameters from the remaining 
parameter space, so that 



P = {Ft\Fi'\p), 

and for every value of p, minimizing 



F{tr;p)^ F^ 



(fc) 



(6) 



(7) 



thereby obtaining a maximum-likelihood estimate of the parameter 
vector p, means that the best-fitting source and background fluxes 
can be expressed in closed analytical form ( lRattenburvll2003h 
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while the non-linear minimization process can be restricted to 
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expressed by means of the magnification rather than the flux, where 

_ „(fc) 

and 

(fe) 

(11) 
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Figure 1. Light curves for a selection of impact parameters mq, wliere 
the event magnification A = (F — Fb)/Fs is plotted as a function of 
{t — t(})/tE, where the closest angular approach uq is realized at to, 
tE = Se/|m|, with 6e denoting the angular Einstein radius, defined in 
Eq. (T), and fi being the relative proper motion between lens and source 
star. 



Figure 2. Light curves for the same impact parameters as adopted for Fig.[T] 
but now the relative flux difference {F — F^,;^^^) / AF with respect to its 
peak is plotted as a function of (t — to)/ti/2, where the half-maximum 
time-scale 4^/2 is defined by Eq. (15). All light curves range between the 
asymptotics for « <S 1 and ii 2> 1 (shown in black), as given by Eq. )19t . 



3 PARAMETERS THAT MATCH OBSERVATIONAL 
FEATURES 

It is well known that avoiding parameter degeneracies and strong 
correlations eases the modelling process. Moreover, a careful study 
of these provides valuable insight into the structure of parame- 
ter space, which can be used as a guide for developing observ- 
ing strategies and explicitly shows their limitations. For ordinary 
microlensing light curves, it turns out that the observed flux can 
be rewritten by means of a different parametrization that better 
matches the observational features, which are key to diagonalizing 
the covariance matrix. 



With the baseline flux i^^aic peak flux Fq''' given by 



F^"^ = F^'^^A[u{to;p)] + Fi'\ 
the maximal flux difference AF''^' reads 



-^0 -^basc 



F^^ [Aiuo) - 1] . 



(13) 



Not only in F*''' and F^^\ but also in F^^^ and AF^'^', the ob- 
served flux F''°' (t; p) is a linear function, namely 



F<'='(t;p) = AF 



(fc) - 1 

A{uo) - 1 



F 



(k) 
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so that these parameters can be separated by a linear fit as well. One 
can then define a half-maximum time ti/2, so that at time to ±ii/2, 
half of the flux offset is encountered, i.e. 



F('^'(to ± ti/2;uo,to,tE)-F^ 



(fc) 



■ AF 



(fc) 
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base 2 

This means that the epochs to ± ^1/2 correspond to a magnification 

Aij2{uo) = A[u{to ± ti/2; U„,to,tE)] = ^ . (16) 

With 111/2 = u{tQ ± ti/2; uo, to, tE,) and by means of Eq. one 
finds 



tl/2 = tE Will 



which leads to 



(17) 
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- 1 



t-tp 
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(18) 



iv in units of the 

base 



If one now plots the offset flux F'''>(t;p) - fHI^ 
difference AF*^*' between peak and baseline, and scales t — to 
with ti/2 rather than tE, as illustrated in Fig.|2] one finds that light 
curves with different uo nearly coincide around their peaks. In con- 
trast, it is the wing region of the event 1.5 ti/2 I'' ~ *o| ^ 3 ti/2 
that is best-suited to provide information about the impact param- 
eter Mo, and thereby of the event time-scale tE, related to the un- 
derlying physical properties^ In principle, mo is determined by the 
slope of the light curve at t = to ± ti/2, but the narrow rang^fl 
strongly limits the feasibility of such an approach in practice. 

It is also instructive to explore the limits of small and large 
separations more closely. Both for it <C 1 and it ^ 1, the observed 
(12) relative flux offset (F'''> (t) - F^^^J/AF**"' becomes independent 
of ito, while ti/2 becomes proportional to ito tE, however with dif- 
ferent proportionality factors for the two extreme cases. As summa- 
rized in Table[T] this is due to the fact that the magnification can rea- 
sonably well be approximated by a more simple expression, namely 
[^(it) — l]/[yl(iio) — 1] approaching uo/u for 11 ^ 1, and (ug/u)^ 
for It ^ 1. For 11 ^ 1, one retrieves the known results that apply 



to microlensing of unresolved sources where F^''' > FJ^'"^ being 
an inevit ability for observations towards M3 1 or other near by other 
galaxies teaillon etalll 19931 : IWozniak & Paczvnski|[T997h . Look- 
ing at Fig. |2] again, one sees that all light curves range between the 
two extreme cases 
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(11 < 1) . 



(19) 



The wings of the light curves converge towards the expression for 



^ The subsequent section shows that the asymptotic degeneracy for large u 
does not extend into this region. 

^ The absolute value of the slope ranges between 0.375 (for uq — * 0) and 
2- 0.586 (for UQ — > CO). 
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Table 1. Asymptotic behaviour of tlie observed flux F^^\ as well as of other relevant quantities, in the limits n <^ 1 or « 2> 1, and its matching parametriza- 
tion. 
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It ^ 1, where for smaller uo, the transition from the w ^ 1 to the 
It 3> 1 asymptotic occurs at larger |t — to 1/^1/2 ■ 



4 DIFFERENT EVENT PHASES AND THE LACK OF 
PREDICTABILITY 

One might think that 3 parameters (like uq, to, and tE) can be ob- 
tained straightforwardly by means of regression from as few as 4 
data points, but such an attempt fails if the observable does not sig- 
nificantly change with a variation of the considered parameter. So 
far, it has been assumed that the light curve has been sampled over 
its full course, so that a full set of characteristics can be determined 
from which model parameters can be derived. However, the light 
curve develops in time, so that some characteristics are not accessi- 

( k) 

ble at early stages. While the baseline flux -F^asc being observed 

(k) 

much before any rise in brightness occurs, the peak flux Fq or the 



flux shift AF('=' = F(S'=' 



^basc' respectively, remain unknown. 



as we shall see more explicitly in the following. 

Reviewing the asymptotics for it <C 1 and u 3> 1 from a 
different perspective already shows us that there is no interpre- 
dictability between the peak and wing regions of the observed light 
curve, and in particular, early observations give poor estimates of 
the peak magnification, as well as on the time-scale tE- In partic- 
ular, maximum-likehihood estimates corresponding to values that 
minimize x^. Eq- freque ntly yie l d very small uo, far away 
from expectations. Therefore. 1 Albrowl j2004h has suggested to use 
a maximum-a-posteriori estimate instead, incorporating the actual 
distribution of the parameters of the observed events of the mi- 
crolensing surveys as prior. 

Let us look into this with a view on the parameter degen- 
eracies. The earliest stages of an event are characterized by |t — 
to\/tE 2> Uo, so that it ~ |i — to|/tE. For the two extreme cases 
of It being far from unity, one finds 
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(ito < It < 1) , 



\t-to\ 

where the second case is realized for a substantial time interval if 
Uo is sufficiently small. 

With the blend ratio p'*^' = F^^''/F^''\ the observed flux 
F'*' {t) can be written as 

"ylKt;p)]+pW" 



= F 
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so that 
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where 
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are characteristic rise times, absorbing tE and g'^''\ with tr = ts — 
tE for g — 0. 

Figure|3]shows F'^''' (t) /f1^2c - 1 as a function of (t - to)/tr 
in a double-logarithmic plot, so that the approximate relations of 
Eq. i22\ correspond to straight lines. It illustrates that different be- 
haviour allows to distinguish three phases of a microlensing event, 
corresponding to the initial rise, a mid-phase, and the peak ap- 
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Figure 3. Development of a microlensing event from a common initial rise 
to different peak fluxes for three selected impact parameters no (colour- 
coded) and two blend ratios g''"' = V^g (solid or dashed hne). 
Double-logarithmic plot showing the relative brightening F''"' (t) / F^'^^ — 
1 above baseline as a function of (t — to )/ tr*"' , where is the rise time 
defined by Eq. )23t , so that for the two selected blend ratios = 
{gW = 0) or tB = 2 tp' (3*'°' = 15), respectively. The black lines 
conespond to the asymptotic behaviour given by Eq. (22). Be aware of the 
fact that the same scale corresponds to tiny changes in the observed flux 
near the bottom of the plot, but huge ones near the top. Similarly, the left 
parts span large time intervals, while the right parts span small ones. 



proach, where a determination of the full model parameter set re- 
quires an assessment of the fundamental characteristics of all these 
phases. With the approximate F'''' (t) diverging as t to, an esti- 
mate for to based on such an approximation corresponds to uq =0. 
Only a departure from the approximation in the form of an evident 
turn-off gives evidence for non- vanishing uq. However, as Fig. [3] 
illustrates, this happens at rather late stage, and roughly when a 
third of the peak magnification has been reached. Given that for 
M < 1, A{u) differs substantially from the asymptotic behaviour 
for u ^ 1, the blend ratio g'*' can in principle be determined 
rather early, which then provides the time-scale ts that is related to 
the underlying physical properties of the event. In practice this re- 
quires sufficiently dense and precise measurements, which are fre- 
quently not available for fluxes close to Jbasc, but the long duration 
of the rise phase in principle allows for lots of data to be collected. 

(k) (k) 

Nevertheless, as soon as both tr and ts can be determined from 
the acquired data, the blend ratio g^''^ and tE are known. Again, 
one sees the power of observations covering the wing of the light 
curve. 

Going the other way round and coming from the peak, the 
microlensing light curve first follows the decay with ti/2, then en- 
ters the mid-phase characterized by ti'"\ and finally follows a de- 
crease described by tP'. Double-logarithmic plots of the relative 
offset brightening ^o'^V^basc — 1 as a function of (t — to)/ti/2 
are shown in Fig. [4] while the relevant magnifications, impact pa- 
rameters, and time-scales for the selected cases are listed in Table[2] 
Again, one sees that the light curves for different blend ratio and tE 
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Figure 4. Development of events from the peak with a given relative bright- 
enmg 

F,5'''Vi^basc' *here light curves corresponding to three selected val- 
ues are shown, towards the baseline flux fP'^ for two different blend 

base 

ratios g = F^^ /F^\ Given that the peak region is characterized by 

its width, ^o^V^basc — 1 shown as function of (t — to)/ti/2 (in ^ 
double-logarithmic plot). As for Figure[3] the asymptotic behavioiff given 
by Eq. j22) is shown by means of black lines. 

can be made to match well near the peak, whereas observations in 
the wing of the light curve on the departure from the asymptotic 
mid-phase behaviour (uo u <^ 1) or on the transition into the 
baseline-approach phase (u S> 1) are useful to determine tE. The 
rather long duration of these phases favours such an attempts by 
not requiring a high sampling rate for obtaining a larger number of 
measurements. 



5 SUMMARY AND FINAL CONCLUSIONS 

Despite the fact that the timescale tE = 6'e/m, where 9e de- 
notes the angular Einstein radius, and ^ the relative proper mo- 
tion between lens and source star, is the crucial one for drawing 
conclusions about the underlying physical properties that led to a 
microlensing event, it is not a good choice for describing the ob- 
servable characteristic features, and thereby not a useful means of 
describing or predicting events in progress. 

Contrary to common belief, 3 model parameters cannot al- 
ways properly be extracted from a least-squares fit involving at least 
4 data points. In fact, such attempts fail if the respective function 
to be matched to the observed data does not significantly depend 
on each of the parameters over the region where data have been 
acquired. Microlensing light curves usually go through 3 phases 
from baseline to peak as well as from peak back to baseline: Two 
rise phases characterized by different rise times tr*"' and ti*' , corre- 
sponding to different power laws of the magnification with the an- 
gular separation between lens and source, as well as a peak region, 
characterized by a half- width ti/2- None of the individual phases 
contains characteristic information about the complete set of model 
parameters, and in order to reveal accurate estimates for all of them, 
each of them requires appropriate coverage. While the mid-phase 
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Table 2. Magnifications, impact parameters and time-scales for the relative peak fluxes adopted in Fig.|4] 



U ' base 




^1/2 


«0 


"1/2 


*E/il/2 


tr /tl/2 


ts /tl/2 


9 = 


50 
20 

5 


50 
20 
5 


25.5 
10.5 

3 


0.0200 
0.0500 
0.203 


0.0392 
0.0956 
0.348 


29.6 
12.3 
3.53 


29.6 
12.3 
3.53 


29.6 
12.3 
3.53 


3 = 15 


50 
20 
5 


785 
305 
65 


393 
153 

33 


0.00127 
0.00328 
0.0154 


0.00255 
0.00654 
0.0303 


453 
177 

38.3 


227 
88.4 
19.1 


28.4 
11.1 
2.39 



tE = 2 4'"' = 16 for g = 15, while = 4*°' = 4*' for 9 = 0. 



gets squeezed for impact angles u 6e of the order of the angular 
Einstein radius 6e or larger, the part of the light curve with u 2> 1 
becomes indistinguishable from the baseline for strongly-blended 
events, so that on the approach to baseline, the offset magnification 
is proportional to u^^ rather than given that still u ^ 1 for 
\t-to\ > 1. 

It is well-known that an appropriate estimate of event pa- 
rameters at early event stages is not feasible, and in particular 
the peak magnification is regularly overpredicted by a maximum- 
likelihood estimate corresponding to minimizing t he sum of nor- 
malized squared deviations x'^- Just for this reason, Albrowl 120041) 
had suggested to use a maximum-a-posteriori estimate instead, with 
a suitable prior. While this brings the estimate closer to its expecta- 
tion value, it does not get around the uncertainty. A closer exima- 
tion shows that the light curve is compatible with an infinite peak 
flux until roughly a third of the true offset magnification is reached. 
Moreover, it is a wing region 1.5 ^ \t ^ to\ ^ 3 ti/2 that 
is best suited to determine the blend ratio ' = F^^ V^s ^^"'^ 
with it the time-scale tE, rather than the immediate vicinity of the 
peak. The rather long duration of this phase allows to obtain a suit- 
able measurement without the need for very dense sampling. 
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For an accurate prediction of the observed flux, a proper deter- 
mination of the full set of model parameters is not required, so that 
local approximations can provide a reasonable substitute. In sharp 
contrast, proper knowledge of tE, which implies knowledge of the 
blend ratio g''^' and the magnification A{t), is a req uirement for 
determ ining the event detection efficiency to planets taaudi et al.l 
l2002h a s well as for prioritising ongoing events in order to maxi- 
mize it ( lHanll2007l . lSnodgrass et al.ff2008l) . Without such informa- 
tion, one neither knows the amplitude, nor the duration, nor the 
location of potentially arising planetary signals. Therefore, an ef- 
ficient campaign for inferring the planet population from observed 
microlensing events needs to invest time into observations that al- 
low to properly determine the event parameters, rather than just try- 
ing to detect planets in poorly determined events, where unsuitable 
assumptions about model parameters may yield to bad choices, or 
efforts could even turn out to be wasted if the planet detection ef- 
ficiency cannot be assessed. Building upon the findings presented 
in this paper, a more detailed study of event (un)predictability tak- 
ing into account the specific capabilities of observing campaigns 
could hence provide important clues towards optimizing strategies 
for detecting planets and determining their population statistics. 



